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Abstract 

We consider a deposition model in which balls rain down at random to- 
wards a 2-dimensional surface, roll downwards over existing adsorbed balls, 
are adsorbed if they reach the surface, and discarded if not. We prove a spa- 
tial law of large numbers and central limit theorem for the ultimate number of 
balls adsorbed onto a large toroidal surface, and also for the number of balls 
adsorbed on the restriction to a large region of an infinite surface. 

1 Introduction and statement of results 

Random sequential adsorption (RSA) is a mathematical model, incorporating stochas- 
tic and geometric elements, for sequential deposition of colloidal particles or pro- 
teins onto a surface (or substrate); particles arrive at random locations, and each 
adsorbed particle occupies a region of the substrate which prevents the adsorption 
of any subsequently arriving particle in an overlapping surface region. Scientific 
interest is considerable; for a series of surveys, see Colloids and Surfaces A, Volume 
165 (2000), for example Privman [12], Senger et al. [14], Talbot et al. [17] and Wang 
[18]. See Evans [5] for a much-cited earlier survey. 

For deposition onto a surface of dimension d, there have been many simulation 
studies, often concerned with the number of particles ultimately adsorbed onto a 
region of substrate. It is of interest to know whether this satisfies a law of large 
numbers (LLN, i.e. a thermodynamic limit) and central limit theorem (CLT, i.e. 
Gaussian fluctuations) as the region becomes large. While previous rigorous math- 
ematical studies were mainly restricted to 1 dimension, for general d Penrose [9] 
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proved a LLN for various continuum systems, and [10] proved both LLNs and CLTs 
for certain lattice systems. Penrose and Yukich [11] proved both LLNs and CLTs 
for continuum systems with finite input where the addition of incoming particles 
is terminated before saturation occurs. Except in the case d — 1 (Dvoretzky and 
Robbins [4]), a CLT for infinite- input continuum RSA remains elusive. 

In the present work we prove a LLN and CLT for an infinite-input continuum 
model related to RSA which has received attention in its own right on grounds of 
realism, namely, a form of monolayer ballistic deposition (BD), representing depo- 
sition in the presence of a gravitational field. Each incoming particle occupies a 
Euclidean ball of radius p in R d+1 , with d = 2 or d = 1; the (d + l)-st coordinate 
represents "height". An incoming particle falls perpendicularly from above towards 
a substrate represented by the surface M d x {0} C R d+1 , which we identify with the 
lower-dimensional space IR d , or a sub-region thereof (the target region). Its down- 
ward motion is vertical until it hits the substrate or one of the particles previously 
adsorbed. If it contacts a previously deposited particle, then the new particle rolls, 
following the path of steepest descent until it reaches a stable position. If the new 
particle reaches the adsorption surface, it is fixed there; otherwise it is removed from 
the system. For d — 1, the model dates back to Solomon ([15], page 129), and the 
formulation in d = 2 by Jullien and Meakin [6] has led to considerable renewed 
interest; see [14, 17]. We state and prove results only for d = 2; changing to d — 1 
makes things easier. 

To avoid having to specify the behaviour of particles near the boundary of the 
target region, we assume, as in most simulation studies, that the target region is 
a torus with integer dimensions. Given A C Z 2 of the form A = {mi, . . . ,m 2 } x 
{m 3 , . . . ,m 4 } (a lattice rectangle), define A C M 2 by 

A = (mi — 1,777.2] x (iTi3 — 1,777.4]. (1-1) 

We focus attention mainly on target regions of this form, and adopt periodic (toroidal) 
boundary conditions for the rolling mechanism. 

Suppose Xi,X 2 , A 3 , . . . are independent and uniformly distributed over A. These 
form the random input to the model with target region A; the vector Aj represents 
the position at which the 7-th incoming ball would end up touching the 2-dimensional 
substrate if it were to fall un-hindered. Successive balls are adsorbed (with possible 
displacement due to rolling) or rejected according to the BD mechanism described 
above, but adopting the toroidal boundary conditions, whereby an adsorbed parti- 
cle near one edge of the rectangle A can influence the rolling of a particle near the 
opposite edge. The process terminates when there is no available space left on the 
substrate large enough to contain a new item (jamming of A); that is, when every 
point of A lies within a distance less than 2p, using the toroidal metric, from some 
point in A that is the location of the point of contact of some previously adsorbed 
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ball. 

Let N(A) denote the (random) number of balls adsorbed at the termination 
time. Our first result is a LLN for N(A) as A becomes large. For any sequence of 
sets (A n ) n >i, set 

liminf(Ai) := U n >i H m > n A m . 

LP 

For p > 1, let — ► denote convergence in p-th moment as n — > oo. 

Theorem 1.1 There is a constant /i = fj,(p) > such that if (A n ) n >i is a sequence 
of lattice rectangles satisfying liminf(A n ) = Z d , then for any pG [1, oo), 

(1.2) 

Next, we give an associated CLT. Let jV(0, a 2 ) denote a normally distributed 
random variable with mean zero and variance a 2 , if a > 0, or a degenerate random 

variable taking the value with probability 1, if a = 0. Let — > denote convergence 
in distribution. 

Theorem 1.2 There is a constant U\ = o"i(p) > such that for any sequence 
(A n ) n >i of lattice rectangles with liminf(A n ) = Z d ; we have as n — > oo that 

\A n \- l V^{N{A n ))^al (1.3) 

and 

\A n \-^(N(A n )-EN(A n )))^ N X0,al). 

Various alternative boundary conditions, other than the toroidal scheme above, 
are also feasible. Particles could simply roll until they touch the surface IR 2 (possibly 
outside the target region); or any particle that ends up touching the surface outside 
the target region could be removed; or (as in Solomon's [15] version of this model, 
generalized to d — 2 in Weiner [19]) the boundary itself could cause a deflection of 
particles (imagine a 'wall' around the boundary of the target region). For the LLN, 
these boundary conditions are not so important, and a result like Theorem 1.1 can be 
obtained any boundary conditions provided the influence of the boundary has finite 
range. Moreover, the target regions in the sequence do not need to be rectangular, 
provided only that they satisfy a condition of vanishing boundary length relative 
to their area. We do not go into details on such generalizations because of their 
proximity to results in [9]. For the CLT, however, alternative boundary conditions 
can cause extra difficulties in the proof. We believe these can be overcome in at 
least some cases of non-toroidal boundary conditions, but have not written out the 
details. 
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While toroidal boundary conditions are usually used in simulation studies, they 
are not so realistic physically. Another way to avoid boundary effects is to take 
the whole of R 2 as target region. Our next result shows that a stationary point 
process, loosely speaking the set of locations of adsorbed points for the BD process 
with target region R 2 , exists as a weak limit of point processes arising from bounded 
target regions. Let S be the space of locally finite subsets of R 2 . For ( e S and 
B C R 2 , let C{B) denote the number of elements of ( in B (so ((■) is a counting 
measure). A point process on IR 2 is a random element ( of S. For more details, see 
for example [3, 16, 13]. 

If ( and ( n (n G IN) are point processes on R 2 , we say the sequence ( n con- 
verges weakly to ( if the finite-dimensional distributions converge, i.e. if for any 
finite collection of bounded Borel sets Bi satisfying ((dBi) = almost surely, the 
joint probability distributions of ( n (Bi) converge weakly to those of ({B^. This is 
equivalent to various other definitions of weak convergence; see e.g. section 9.1 of 
[3]. 

Given a lattice rectangle A, let £ A be the point process of locations in R 2 of 
ultimately accepted particles, for the BD model with target region A (which will be 
a point process in R 2 , all of whose points lie in A). Our next result concerns weak 
convergence of the point process £ A as A becomes large. As with Theorem 1.1, the 
result is not sensitive to the toroidal boundary conditions. 

Theorem 1.3 There exists a stationary point process such that if (A n ) n >i is 
any sequence of lattice rectangles with lim inf n ^oo A n = Z 2 7 the sequence of point 
processes £ An converges weakly to £. 

Given any bounded region B C R 2 , the interpretation of £,{B) is as follows. The 
variable (, A (B) is the number of adsorbed points in B when the target region is A. 
As A becomes large this has a weak limit which is ^(B). If one now, in turn, makes 
B large, it is of interest to know if C,(B) satisfies a CLT, and our final result says 
that this is indeed the case. We restrict attention to rectangular regions although 
other shaped regions can also be dealt with (see the proof). 

Theorem 1.4 There exists a constant a 2 = o"2(p) > such that for any sequence 
(_£>„)„>! of lattice rectangles with lim inf (.£?„) = Z d , we have as n — > oo that 

l^l-Var^))^ 2 (1.4) 

and 

\B n \- 1/2 [t(B n ) - Et(B n j) A/-(0,a 2 ). 
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Other properties of the point process £ are also of interest. The proof of Theorem 1.4 
involves showing that £ has exponentially decaying correlations, which is of interest 
in its own right. 

Weiner [20] considered an alternative version of the BD model in two (or more) 
dimensions in which the region of substrate occupied by a particle is a rectilinear 
square rather than a circle (the 'Solomon model'). He claimed to prove a CLT for 
the number of particles ultimately adsorbed onto a large target region. However, his 
argument uses assertions from Weiner [19], which he later retracted (Weiner [21]). 
It is possible to adapt our methods to yield a CLT for Weiner 's 'Solomon model', 
at least in the case of toroidal boundary conditions, partially vindicating Weiner's 
claims regarding this model (though not the 'Renyi model'). 

2 Geometric preliminaries 

In notation used throughout this paper, denotes the origin (0, 0) of M 2 . For 
x = (xi,X2) G IR 2 , \\x\\ denotes the Euclidean norm (modulus) \fx\ + x\ of x. For 
x G 1R 2 and R C M 2 , x + R denotes the translated set {x + y : y G R}. For r > 0, 
define the continuum disk and D(x; r) C M 2 the lattice ball B(x; r) C Z 2 by 

D(x; r) — {y : \\y — x\\ < r}; B(x; r) = D(x; r) fl Z 2 . 

If E is an event in a given probability space let 1e be the indicator random variable 
taking the value 1 if E occurs and if not. Finally, for any directed graph, by a root 
of the directed graph we mean a vertex with indegree zero. 

We start with two purely geometric results about the mechanics of the BD model 
with target region given by the infinite surface M 2 x {0}. Each accepted particle lies 
on the substrate, and so can be represented by the point in M 2 at which it touches 
the surface. The position of an accepted particle is a translate (or displacement) of 
the location in M 2 above which it originally comes in. The displacement (and also 
the decision on whether or not to accept) is determined by the initial location at 
which the particle comes in, and the positions (after displacement) of the previously 
accepted particles. 

Lemma 2.1 With probability 1, no particle receives a displacement of modulus 
greater than 8p. 

Proof. Choi et al. [2] enumerate the possible fates an incoming ball might 
undergo. Since these involve at most 4 deflections, in effect they state the result but 
do not give a complete proof. Therefore we do so here. 
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Let each particle already accepted be represented by a point in M 2 located at its 
point of contact with the substrate. For any two such points the inter-point distance 
r (say) satisfies r > 2 p. 

For a new particle, let it too be represented by a point in M 2 , obtained by project- 
ing the position of its center down onto the substrate (imagine looking down on the 
substrate from above). As it rolls, the point representing the new particle performs 
a piecewise linear motion in M 2 . The first line segment of this motion represents 
an initial period when the new particle touches a single existing particle, and is of 
length at most 2p. Thereafter, each successive line segment will represent motion 
while touching two existing particles, and will be along the mediator (perpendicular 
bisector) between the two points representing those particles. Let dj denote the 
distance between the two particles which the new particle touches during the j-th 
step of linear motion, and note that 2p < dj < 4p. 

Each change in direction of this piecewise linear motion in M 2 , say from step j to 
step j + 1 of linear motion, will occur at the circumcenter of three existing points. 
If this circumcenter lies inside the triangle with corners at those three points, then 
the motion comes to a stop and the particle is discarded, according to the BD rules. 
Therefore for the motion to continue, the circumcenter lies outside this triangle, so 
that the triangle must have an obtuse angle. The inter-point distance dj + \ is the 
longest edge length of this triangle, while dj is one of the other two edge lengths. 
Since the triangle has an obtuse angle, and all three edges are of length at least 2p, 
we obtain 

d) +l >d] + Ap\ 

and since the edge-lengths dj must all be at most 4p, this means that the sequence 
(dj) terminates in at most three steps, in addition to the initial rolling in contact 
with just a single previous particle. Since each piecewise linear step is of length at 
most 2p, this completes the proof. | 

Let us say that after k adsorptions, a given point x G M 2 is available for a particle 
to be adsorbed, if there are no adsorbed particles touching the surface in D(x; 2p). 

Lemma 2.2 There exist e > 0, £i > with the following property. Suppose that 
for some k, after k adsorptions, a given point x G M 2 is available for a particle to 
be adsorbed. Then there exists a region of area at least e± such that any incoming 
particle with location in that region will be adsorbed in a position that makes all 
points in D(x;eq) unavailable. 

Proof. Assume without loss of generality that p = 1/2 and x = 0. Take sq < 1/8. 
First suppose no adsorbed particle lies within distance 1 + e of the origin. Then 
any particle arriving in the ball D(0;e ) will be accepted without rolling, and for 
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such a new particle the unit diameter ball centred at that particle covers the ball 
D(0;e ). 

Next suppose there already exists a particle (at x, say) with 1 < ||x|| < 1 + Eq. 
A particle arriving in D(0;eq) fl D(x; 1) will receive a first deflection and roll, but 
not very far. This is because its initial distance from x is more than 1 — eq so it 
initially rolls at most a distance e before it reaches the surface or receives a second 
deflection. If a second deflection takes place, at that instant the new particle then 
lies on the mediator of two adsorbed points x, y say. The distance between x and y 
is more than 2(1 — e ). If a third deflection were to take place it would be at the 
circumcenter of adsorbed points points x, y, z, say, making a triangle with an obtuse 
angle. But this is impossible; for example if (y, z) were the longest edge, the cosine 
rule would give us 

\\y - z\\ 2 > \\x - z\\ 2 + \\x - y\\ 2 > 1 + (4 - 8e + 4^), 

and therefore \\y — z\\ > 2 and the third deflection does not take place. Therefore 
after the second deflection the linear motion terminates either with adsorption or 
rejection. By Pythagoras' theorem, the distance travelled in this last linear motion 
after the second deflection is at most 

V / l-(l-^o) 2 <(25 ) 1/2 

and therefore the total displacement of the particle before adsorption is at most 
5 + (2e ) 1/2 . 

Therefore if a particle arrives within a distance Eq of the origin, it is adsorbed or 
rejected at a distance at most 2e + (2£o) 1 ^ 2 from the origin. Since Eq < 1/8, this 
is at most 3/4. If adsorbed, it therefore prevents any subsequent adsorption taking 
place in D(0; e ). Therefore we are done, unless there is a possibility of rejection for 
particles arriving in D(0;e ). 

Next, suppose that it is possible for a particle arriving within distance £o of the 
origin to be rejected. If this happens it will be at the circumcenter of points x, y, z 
(say) after initial deflection by x and subsequent deflection by y (say). In this case 
the circumcenter of x, y, z is at a distance less than 1 from each of x, y, z, and every 
point inside the triangle xyz is unavailable. In particular, the origin does not lie 
inside the triangle xyz; however, it does lie within distance 2e of the midpoint of x 
and y. 

Now suppose a particle lands in D(0; e ) on the same side of the line xy as the 
origin. In the course of its subsequent rolling it stays on the same side of the line 
xy as the origin; if not there would be some other line segment between adsorbed 
centres, other than {x, y}, of length between 2(1 — e ) an d 2, whose midpoint lies 
in -D(0; 2e ), which is geometrically impossible provided £ is sufficiently small. At 
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the end of its motion, if it were rejected, that would take place at the circumcenter 
of points xyz', say, and in that case all points landing in xyz' would be unavailable. 
However, provided Eq is small enough, the origin must lie in xyz', and therefore we 
would have a contradiction. 

It follows that provided e is small enough, a point landing in D(0;eo) on the 
same side of the line xy as the origin will be accepted, in a position that makes the 
region D(0;e ) unavailable. The desired conclusion follows, with e\ = tt£q/2. | 

3 Probabilistic preliminaries 

The author [8, 10] has developed general LLNs and CLTs for functionals on the 
restriction of spatial white noise processes to finite regions of the lattice, as follows. 
Suppose (E,£,P ) is an arbitrary probability space, and X = (X x ,x G Z 2 ) is a 
family of independent identically distributed random elements of E, each X x having 
distribution P . Let X' be the process X with the value X at the origin replaced by 
an independent copy X* of X (that is, an _E-valued variable X* with distribution 
P , independent of X), but with the values at all other sites the same. Let TZ 
denote some collection of nonempty finite subsets of Z 2 , with x + B G TZ for all 
B e K,x e Z 2 . 

A stationary TZ-indexed functional of X is a family H = (H(X; A), A G TZ) of 
real-valued random variables, with the property that (X x ,x G A) determines the 
value of H(X; A) and H(r y X; y + A) = H(X; A) (almost surely) for all y G Z 2 , 
where r y X is the family of variables (X x - y ,x G Z 2 ). Let A (A) be the increment 
H(X; A) — H(X'; A). The functional H stabilizes on sequences tending to Z 2 if there 
exists a random variable A (oo) such that for any 1Z- valued sequence (Ai)n>i with 
liminf n ^ 00 (A n ) = Z d , the variables A (A n ) converge in probability to A (oo). 

A stationary TZ-indexed summand is a collection (Y Z (X;A),A G 7Z, z G A) of 
real-valued random variables with the property that (X x ,x G A) determines the 
value of Y Z (X;A), and Y y+z (T y X;y + A) = Y Z (X; A) (almost surely) for all y G 
Z d ,A G 7Z,z G A. The associated induced stationary TZ-indexed sum is given by 
H(X; A) = ^2 A Y Z (X; A), which is a stationary 7£-indexed functional. 

We restrict attention here to the case where TZ is the collection, here denoted 
£>, of all lattice rectangles {mi, . . . , m 2 } x {m 3 , . . . , m 4 } with m 2 — m x > 20p and 
m 4 — m 3 > 20p. This is different from the class of sets denoted B in [10]; nevertheless, 
the following general law of large numbers is proved in just the same manner as the 
first part of Theorem 3.1 of [10]. 

Lemma 3.1 Suppose (Y Z (X; A) : A G B, z G A) is a stationary B-indexed summand 
inducing a stationary B-indexed sum (H(X; A); A G B). Suppose that 

sup{E|Yo(X; A)\ : A G B : G A} < oo. (3.1) 
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Suppose there exists an integrable random variable Y (X) such that Y (X;B n ) — > 
Y (X) in L 1 as n — > oo ; for any B-valued sequence (B n ) n >i with liminf(_B n ) = Z 2 . 
If (A n ) n >i is a B-valued sequence with liminf(A n ) = Z 2 , then 

\A n \~ 1 H(X; A n ) EY as n -> oo. (3.2) 

Let JF be the cx-field generated by (X y ,y ^ 0), where y ^ means ?/ precedes 
or equals in the lexicographic ordering on Z 2 . The following general CLT is a 
corollary of Theorem 2.1 of [8] (see Remark (iii) thereafter in [8]). 

Lemma 3.2 Suppose (H(X; A); A G B) is a stationary B -indexed functional of X 
which stabilizes on sequences tending to Z 2 , and for some 7 > 2 satisfies 

supE[|A (A)n < 00. (3.3) 

AeB 

Suppose that (Ai)n>i is a B-valued sequence with liminf(A„) = U 1 . Then as n — > oo ; 
\A n \~ l Xax(H(X; A n )) converges to a 2 , and 

lA^iH^An) -EH(X;A n )) -^^(O,. 2 ), (3.4) 

with a 2 =E[(E[A (oo)|Jo]) 2 ] . 

We also use the following CLT for stationary random fields, from Bolthausen [1] . 
For A U A 2 C Z 2 , let d(A 1 ,A 2 ) = inf{(||^i - z 2 \\ : z { G A u i = 1,2}. Let dA x be the 
set of z G Z 2 \ A 1 such that d(A 1 , {z}) = 1. 

Lemma 3.3 [1] Suppose (ip x ,x G Z 2 ) is a real-valued stationary random field. For 
integers ai,a 2 ,n > 1 define 

) = sup{|P[Fi fl F 2 ] — P[Fi\P[F 2 ]\ : F; G a{jp z : z G Aj), 

|A| < a i: d(A u A 2 ) > n}. 

Suppose Y^m=i ma ai,a2( m ) < 00 / or a i + °2 < 4, and ai i00 ( m ) = o{m~ 2 ), and 
EW^o?] < oo, and Zlm=i ma i,i( m ) 1/3 < 00 ■ 

Then a 2 := XLez 2 Cov(-0o, ^z) converges absolutely, and if a 2 > 0, then for any 

sequence (T n ) n >i of subsets of 7? with\d(T n )\/\T n \ ^ asn^ oo, |r n |~ 1 / 2 XLer 
A'(0.n 2 ). 
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4 Proof of LLN 



Let V be a homogeneous Poisson point process of unit intensity on M 2 x [0, oo). Given 
AeB, label the points of the restriction of V to A x [0, oo) as {(X(A), T^A))}^ 
with Ti(A) < T 2 (A) < T 3 (A) < Throughout the proofs of Theorems 1.1-1.3, 
we assume without loss of generality that the random input for the variable N(A), 
defined in the introduction, is given by the sequence of variables Xi(A), X 2 (A), . . . 
representing the locations of successive incoming particles (thus Ti(A) is taken to be 
the time of arrival of the i-th incoming particle). By this device, coupled realizations 
of N(A) are defined for all A G B simultaneously. 

For each point (X, T) of the restriction of V to A x [0, oo), define the pair 
/(X, T; A) = (J (X, T; A), I^(X, T; A)), with I (X, T; A) G {0, 1} and I^(X, T; A) G 
M 2 , as follows. Let Iq(X,T;A) (an indicator variable) be equal to 1 if the ball ar- 
riving at location X at time T is accepted, and to zero if it is rejected, in the 
realisation of the BD model with target set A described above. If I (X,T; A) — 1, 
let 7_>(X, T; A) denote the lateral displacement received by the particle arriving at 
X at time T, prior to being adsorbed. If Io(X, T; A) = 0, let i_+(X, T; A) = 0. 

By Lemma 2.1, the decision on whether to accept an incoming particle, and also 
its displacement if accepted, are determined by the locations (after displacement) of 
those particles lying within a distance at most lOp from the location at which the 
new ball arrives. 

The proof of the LLN and CLT involves a graphical construction. Make the 
Poisson process V on IR 2 x [0, oo) into the vertex set of an (infinite) oriented graph, 
denoted Q, by putting in an oriented edge (X, T) — > (X',T') whenever \\X' — X\\ < 
20p and T < V . Observe that particle (X, T) cannot affect (X', V) directly unless 
there is an edge (X, T) — > (X', T') or the toroidal boundary conditions come into 
play. 

For z — (zi, z 2 ) G Z 2 , and e > 0, define the squares 

Qz,e ■= {{zi - l)e, zis] x ((z 2 - l)e, z 2 e]; Q z := Q zA . 

For x,y G Z 2 , let us say that y is affected by x before time t if there exists a 
(directed) path in the oriented graph that starts at some Poisson point (X, T) with 
X G Q x , and ends at some Poisson point (Y,U) with Y G Q y and U < t. Let 
E t (x,y) denote the event that y is affected by x before time t. 

Lemma 4.1 There is a constant Si G (0, 1) such that for all x, y G 7L d , 

P[E Slllx - y]l (x,y)]<2(3-\\*-vH). 

Proof. See Lemma 3.1 of [9]. This applies directly if 20 p < 1, and its proof is easily 
adapted to the case 20p > 1. | 
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For z <EZ d and t > 0, define the cluster C z>t C Z d by 

Cz,i := G Z d : 2 is affected by x before time t}, (4.1) 

which is almost surely finite by Lemma 4.1 and the Borel-Cantelli lemma. Let B zt 
be the smallest element of B that contains U xg c* z t B(x; 4 + 20p). Note that 5 Z; t 
includes a 'buffer zone' around C z j so that 

dist {u ye i.d\ Bzit Qy, U xeCz , t Qx) > 20p, 

so that even if we were to add extra points outside the union of squares Q y , y G B Zjt , 
there will not be any connected path in the graph from any of these added points 
to any Poisson point (X, T) £ Q z x (0, £]. This will be important later on. 

Lemma 4.2 Suppose z G l?,t > 0. If A is a lattice box with B zt C A, £/ien for all 
Poisson points (X, T) lying in Q z x [0,t] we have I(X,T; A) = I(X,T; B Zyt ). 

Proof. By definition, the influence of Poisson points outside B z t does not propagate 
to any Poisson points in Q z x [0,t]. Therefore the fate of such points is the same 
whether the target region is A or B z>t . g 

Let So be the space of all finite subsets S of D(0; lOp) such that \\x — y\\ > 2p 
for all distinct x,y G S. Define : —> {0, 1} and : So — > IR 2 as follows. For 
S G So, let \E f o(5') take the value 1 (respectively 0) if an incoming particle at the 
origin of M 2 is accepted (respectively rejected), given that S is the configuration of 
previously accepted particles in B(0; lOp). If V (S) = 1, let G M 2 be the 

(lateral) displacement of an incoming particle at the origin of M 2 , prior to acceptance, 
given that S is the configuration of previously accepted particles in -5(0; lOp). If 
*o(S) = 0, let *^(S) = G R 2 . 

We construct a spatially homogeneous form of the BD model with the whole of 
M 2 as target region, as follows. Define subsets Gi, i = 0, 1, 2, 3, . . . , of V as follows. 
Let Go be the set of roots of Q, and recursively, if Go, • • • , G ra are defined, let G n+ i 
be the set of roots of the graph Q with all vertices in Go, • • • , G n removed. As a 
consequence of Lemma 4.1, the sets G ,Gi,G 2 , . . . form a partition of V (see [9], 
Lemma 3.2). 

Define I(X, T) = (J (X, T), J_(X, T)) with I (X, T) G {0, 1} and T) G M 2 

(representing acceptance status and lateral displacement respectively) for (X, T) G 
G ,G 1? ... as follows. If (X,T) G G then set I (X,T) = 1 and I^(X,T) = 0. 
Recursively for n — 1, 2, 3, . . . , for (X, T) G G n , set 

S X , T = {Y + I^(Y, U)-X: (Y, [/) G U^G,,,, J (F, C/) = 1, 

||y + /_(y,c/) -x|| < i0p}, 



n 



then set I (X,T) = V (S x ,t) and I^(X,T) = ^^(S x ,t)- 

For t > 0, let £ t be the point process of positions after displacement of particles 
accepted up to time t; that is, re-labelling the points of V in arbitrary order as 
{(XjjTj)}^!, let £ t be the random locally finite set in M 2 defined by 

6 = {X 3 + L4X„ Tj) : I (X j: Tj) = 1, T, < f}. (4.2) 

This point process is now rigorously defined in terms of the Poisson process V and 
the graphical construction. It is a stationary point process on M 2 . Define the limiting 
point process 

i = U t >„&. (4.3) 

Similarly, for A G B let £ 4 A be the set of locations after displacement (using the 
toroidal boundary conditions) of points (X, T) of V fl (A x [0, oo)) such that T <t 
and /o(^) ^4) — 1- All points of £f lie in A. Define the limiting point process 
£ A = U t >o£ t A , a point process in A. 

Choose 62 G (0, Sq/2), with £o taken from Lemma 2.2, and with I/E2 G IN. For 
z G let /3 2 denote the (random) time at which the square Q Zy£2 becomes blocked, 
i.e. the first time at which the point process leaves no part of the surface of the 
square Q Zt£2 still available to adsorb a sphere. For A G £>, define j3 A in the same 
way with respect to the point process £f. The next result says that the variables 
j3 z and (3£ are almost surely finite and in fact their distributions have tails which 
decay exponentially, uniformly in z, A. 

Lemma 4.3 It is the case that 

lim sup t' 1 sup {log P\j3 g > t] } < 0, (4.4) 

and 

limsupr^upflogP^ >t]:AeB,ze Z 2 ,Q Z , £2 C A} < 0. (4.5) 

Proof. Suppose f3 z > t, i.e. at time t there exists a point x G Q z ,e 2 that is not 
yet covered, i.e. still available. Then by Lemma 2.2, the probability of a particle 
arriving and causing Q Zj£2 to be covered by time t+h is at least eih+o(h). Therefore 
we can choose hi > such that 

P[f3 z >t + h 1 \/3 g >t] < 1- £i/ii/2, 

so that by induction, P[(3 Z > nhi] < (1 — Eihi/2) n for all n G INT. This argument 
holds uniformly in 2, and (4.4) follows. Furthermore, the same argument carries 
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through to the torus, to give (4.5). | 
For y G Z d , let 

J y = max{/3 2 : Q z>£2 cQ x ,xE B(y; 4 + 20p)}, (4.6) 

and (for y e A e B, with || • ||a denoting the toroidal metric) 

J* = max{/3? : Q z , £2 C Q x ,x e A,\\x - y\\ A < A + 20p}. (4.7) 

For each x G B(y; 4 + 20p), the square Q x is jammed by the point process £j y , 
meaning that it is not possible for any Poisson point arriving after time J y to be 
accepted at a position in Q x . In particular, by Lemma 2.1, all particles arriving in 
Q y after time J y are rejected. Define the enlarged 'cluster' B' y by 

B' y = U xe B(y;4+20p)B x ,J y . 

Using this enlarged cluster we can strengthen Lemma 4.2 to account for arrivals at 
all times, as follows. 

Lemma 4.4 Suppose y G Z 2 . If A is a lattice rectangle with B' y C A, then for 
all Poisson points (X,T) lying in Q y x [0, oo) we have I(X,T;A) = I(X,T; B' ) = 
I(X.T). 

Proof. Suppose (X, T) is a Poisson point in Q x x [0, J y ], for some x G B(y; 4 + 20p). 
Then by Lemma 4.2, we have I(X, T; A) = I(X, T; B' y ) = I(X, T). 

By definition of J y , it follows that the restriction of the point process £j to 
the set U X £B(y;2+w P )Qx precludes the subsequent adsorption of any more particles in 
^xeB(y;2+w P )Qx, and in particular prevents acceptance of any subsequent particles 
arriving in Q y ; therefore for every Poisson point (X, T) G Q y x (J y ,oo), we have 
I(X, T; A) = I(X, T; B' y ) = (0, 0). | 

By Lemma 4.3, the variable J y has an exponentially decaying tail, uniformly in 
y, i.e. 

limsupt -1 sup log P[J y > t] < 0. (4.8) 

For z G Z d , let X z be the image of the restriction of V to Q z x [0, oo), under the 
translation that sends each point (X, T) to (X — z, t). This is a homogeneous Poisson 
point process on Q x [0, oo). The Poisson processes {X z ,z G Z 2 ) are independent 
identically distributed random elements of a measurable space (E,S), where E is 
the space of locally finite subsets of Qo x [0, oo). The idea behind the proof of 
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Theorems 1.1 and 1.2 is to regard N(A) as a stationary £>-indexed functional driven 
by the process X = (X z ) zeZ 2, and use Lemmas 3.1 and 3.2 from Section 3. We write 
X rather than X in this case. 



Proof of Theorem 1.1. For A e B and z G Z 2 , 



set 



Y Z (X;A) 



E 



J(X,T;A). 



(4.9) 



(X,T)£Vn{Q z x{0,oo)) 



Then (Y Z (X; A), A <E B, z <E A) defined by (4.9) is a stationary £>-indexed summand 
on on the process X = (X z ) zeZ 2, and the corresponding stationary £>-indexed sum 
H(X; A) is equal to N(A). It suffices to check the conditions in the general result 
Lemma 3.1. Since the variables Y Z (X; A) are uniformly bounded by a constant, (3.1) 
holds. 

Suppose (Aj) n >i is a B- valued sequence with liminf(A n ) = Z d . Then there ex- 
ists a random variable N\ such that for all n > JVi, B' C A n . By Lemma 4.4, for 
all n > N\ and all Poisson points in Qo x [0, oo) we have /(X, T; A n ) = /(X, T; 5q). 
Hence Y (X; A n ) = Y (X; B' ) for n > Ni, so Y (X; A n ) converges almost surely to 
Y (X;B' ) as n — > oo. Therefore Lemma 3.1 applies to give us (1.2) with conver- 
gence in L 1 . Convergence in LP then follows since N(A)/\A\ is uniformly bounded 
by a constant, g 

Proof of Theorem 1.3. Suppose (A n ) n >i is a i3- valued sequence with liminf(yl n ) = 
Z 2 . It suffices to prove that for any bounded Borel B C W 1 we have almost sure 
convergence 



If n is sufficiently large so that B' z C A ra for all z within distance 4 + lOp of 5, 
then by Lemma 4.4, £ An (B) = £(B), which gives us (4.10). g 



With as defined in the previous section, (N(B), B e i3) is a stationary i3-indexed 
functional driven by the white noise process A\ Our goal is to apply Lemma 3.2. 

In this setting, the process X' appearing in the conditions for Lemma 3.2 is 
obtained from the process V 1 defined as follows. Let X' be the Poisson process 
obtained by replacing the restriction X of V to Q x [0, oo) with an independent 
Poisson process X* on Q x [0, oo), so that 



^»(£)-£(£J). 



(4.10) 



5 Proof of CLTs 



X' = (V\(Q x [0,oo)))U^. 



(5.1) 
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The points of V'\(Qo x [0, oo)) are the same as those of V\(Qo x [0, oo)). However, 
the decision on whether to accept may be different; let I'(X, T) = (I' (X, T), IL>(X, T)) 
be defined in the same manner as I(X, T) but based on the process generated by V' 
rather than V; likewise, given A G B, for (X, T) G V'D (A x [0, oo)) let I'(X, T; A) = 
(r (X,T; A), I'_^(X,T; A)) be defined in the same manner as I(X,T;A) but based 
on the process generated by V rather than V. 

Lemma 5.1 Let z G Z d , and t > 0. Suppose £ B' z . Then I(X,T;A) = 
I'(X, T; A) for all points (X, T) ofV in Q z x [0, oo), and all A e B with B' z C A. 

Proof. Suppose (X,T) is a point of V, with X e Q z . Then 

/(X, T; A) = I(X, T; B' x ) = I'(X, T; B' z ) = I'(X, T; A), 

where the first equality comes from Lemma 4.4, and the second comes from the 
equality of V and V outside Qo x [0, oo). | 

The idea for proving stabilization goes as follows. By Lemma 4.1, the effect of 
changing the inputs at the origin propagates like an 'infection' spreading through 
space at a linear rate. However, if this 'infection' encounters a 'wall' of thickness 
lOp surrounding the origin, consisting of sites which are entirely blocked before the 
infection reaches them, then this wall prevents any spread of the effect of changing 
inputs at the origin to the other side of the wall. The existence with high probability 
of a wall surrounding the origin follows from the fact that the probability that a site 
is not blocked by time t decays exponentially in t. 

By Lemma 4.1 and (4.8), P[0 G B' z \ decays exponentially in ||z||. Therefore, by 
the Borel-Cantelli lemma, 

P[0 G B' z for infinitely many z] = 0. (5.2) 

For t > 0, define the annular region 

•At '■= U zeB{O;t+A+20p)\B{Q;t)Q z 

and the 'distant' set 

T^t '■= ^zeZ 2 \B(O;t+A+20p)Qz- 

Lemma 5.2 Let t > 0, and suppose £ B' z for all z G A t n Z d . Then I(X, T; A) = 
I'(X, T; A) for all points (X, T) ofV in A t U V t , and all A G B with B' z C A for all 

zeA t nz d . 
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Proof. Suppose ^ B' z for all z G At- Suppose also AeB with B' z C A for all 
z G At- By Lemma 5.1, 

I(X,T;A) = I'(X,T;A), V(X, T) G P n (A t x [0, oo)). (5.3) 

Next consider Poisson points (X, T) with X G T> t . Any path to (X, T) from 
Qo x [0, oo) must pass through At x [0, oo), and the status of all points in this region 
is unaffected by the change in Vq, so that /(X, T; A) = I'(X, T; A). More formally, 
we use an induction, as follows. 

Define generations G (A, i), G\{A, t), . . . as follows. Let Q(A, t) be the restriction 
of Q to vertices in (V t C\A) x [0, oo). Let Go(A,t) be the set of roots of Q(A,t). Then 
remove vertices of Go(A,t) from Q(A, t) and let Gi(A,t) be the set of roots of the 
remaining oriented graph. Then remove vertices in G\(A,t) too and let G2{A,t) be 
the set of roots of the remaining graph and so on. The sets G (A, t), Gi(A, t), . . . 
form a partition of the vertex set of G(A, t), because C z j defined at (4.1) is finite 
for all z, t. 

The inductive hypothesis is that if (X, T) G G n (A, t), then /(X, T; A) = I'(X, T; A). 
This is true for n — 0, because if (X, T) G G (A, t), then any (X', T') for which X' G 
A and there is an edge from (X',T') to (X, T) lies in the annulus A t x [0, oo), and 
therefore by (5.3) satisfies I'{X',T';A) = I(X',T';A), which implies /'(X, T; A) = 
I(X,T; A), since the decision on the value of J(X, T; A) depends only on the deci- 
sions at points (X', T") from which there are edges to (X, T) 

Now suppose the hypothesis is true for n — 0, 1, . . . , k — 1. Then if (X, T) G 
Gk(A,t) all of the points (X',T') from which there is an edge to (X, T) lie either 
in one of the generations G (A,t), . . . ,Gk~i(A,t), or in At x [0, oo), and there- 
fore, by the inductive hypothesis and by (5.3), all such (z',T') satisfy I(z',T'; A) = 
I'(z', T'\ A) , and hence again I'(z, T; A) = I(z, T; A). | 

Lemma 5.3 If for AeB we set 

A o(^) = E E (I(X,T;A)-I'(X,T;A)) 

zeA (x,T)er-.xeQz 

then 

supE[|A (A)| 3 ] < oo. (5.4) 

AeB 

Proof. Modify the graphical construction from Section 4 to produce an oriented 
graph Q A with vertex set the Poisson process V H (A x [0, oo)), by putting in an 
oriented edge (X, T) — > (X',T") whenever ||X' — X||^ < 20p and T < T", where 
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|| ■ || a denotes the toroidal metric. For x, y G A, let us say that y is affected in A by 
x before time t if there exists a (directed) path in the oriented graph Qa that starts 
at some Poisson point (X, T) with X G Q x , and ends at some Poisson point (Y, U) 
with Y G Q y and U < t. For z G A and t > 0, define the 'A-cluster' C^ 4 to consist 
of all x G A such that some y £ A with ||y — z\\a < 4 + 20p is affected in A by a; 
before time J^. 



A similar argument to the proof of Lemma 5.1 yields I(X,T;A) = I'(X,T; A) 
for all (X, T) with X G Q z and 5(0; 4 + 20p) n = 0. Hence for any A E B, 



However, P[_B(0;4 + 20p) fl ^ 0] decays exponentially in ||z||, uniformly over 
A G B, because of Lemmas 4.1 and 4.3. The bounded moments condition (5.4) 
follows. I 

Proof of Theorem 1.2. As in the proof of Theorem 1.1, define Y Z (X\ A) by (4.9). 
The aim is to apply Lemma 3.2. The bounded moments condition (3.3) follows from 
Lemma 5.3. 

We need to check stabilization. By (5.2), there exists a random R such that 
^ B' z for z G Z d with ||z|| > R. Let -B^ be the smallest element of B containing 
^zeB(o-,z)B' z . Suppose (Ai)n>i is a i3-valued sequence with liminf(A n ) = Z d . Then 
there exists random N 2 such that 



and such that a similar expression holds with regard to the Poisson process V' rather 
than V . Then by Lemma 4.4, for all n > N 2 and for all z G B(0; R), we have 



Y Z {X- A n ) - Y Z {X'- A n ) = Y Z {X- B^) - Y Z {X'- B^) for all n > N 2 . (5.5) 




Boo C A^, for all n > N 2 . 



Also, there exists random N 3 > N 2 such that for all n > N 3 we have 



5; C A n for all 2 G Ar fl Z d . 
Hence by the definition of R and Lemma 5.2, for all z G Z d \ -8(0; -R), 



A,) = F 2 (A"; A n ) for all n > N 3 . 



Combined with (5.5), this gives us for all n > N-. 



(Y Z (X; A n ) - Y Z (X>; A n )) = ^ (Y Z (X ; B^) - Y Z (X'; B^)), 
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which demonstrates stabilization of the induced functional 

H(X;A) = Y t Y v (X;A) = N(A). 

yeA 

Therefore all the conditions for Lemma 3.2 hold here, and by that result the con- 
clusion of Theorem 1.2 holds, subject to showing that o\ > 0. 

The value of <7\ is independent of the choice of sequence (A n ) (provided lim inf (A n ) = 
Z 2 ) and therefore to show <ji > using (1.3) we are at liberty to choose any sequence 
(A)n>i- Let K = |~200p~|. Take A n to be a lattice square of side Kn, and divide 
A n into squares of side K, which we shall refer to as blocks. 

Inside each block Si let % be the annulus of thickness 18p, consisting of points 
at a distance more than 2p but less than 20p from the boundary of the block. Also 
let S~ be the interior region consisting of points at a distance more than 20p from 
the boundary of the block. Let Ii be the indicator random variable of the event that 
before there are any arrivals at all in Si \ Ti , a sequence of ball centres arrive in in 
such a way that the corresponding particles are adsorbed without rolling and cause 
a barrier between the interior region S~ and the complement of Si, by making all 
points in Tj unavailable. 

The probability P[ij = 1] is very small but not zero, and does not depend on i 

2 

or n since K is fixed. Let N n = X^=i ^- Then E [iV n ]/| A| is a non-zero constant. 

Let T be the cr-field generated by the value of N n , along with the positions of 
the accepted particles not lying in the union of the squares {S~ : Ii — 1}. Then 

Var(if (AO) = Var(E [H{A n )\F\) + E [Var(#(A|^))] > E [Var(//(A|^))]- 

Suppose we are given the value of N n and the configuration of accepted items out- 
side the squares S~ , Ii = 1. The only remaining variability is from the number of 
accepted items inside the inner squares S~ contributing to N n . 

Let S* be the square consisting of points in Si at a distance at least 22p from 
Si, i.e. a slightly smaller square inside S~. We consider two possible lattice config- 
urations inside S*. 

Let {xi t i, . . . , Xi >ni } be the restriction to S* of a regular triangular lattice with 
each point distant (2.02)p from its neighbours (tight packing). Let {yi,!, . . . , 2/i, n2 } 
be the restriction to S* of a regular triangular lattice with each point distant 3p 
from its neighbours (loose packing). Let Ei be the event that the first ri\ particles 
in Sj~ to arrive are in the small disks D(xi, 0.01), 1 < % < ri\, and are in different 
disks. Let Fj be the event that the first n 2 particles in S~ to arrive are in the disks 
D(yi, 0.01), 1 < % < ri2, and are in distinct disks. 

Events Ei and F$ have probability bounded away from zero. More particles 
will be packed into the square S~ on event Ei than on event Fj. It follows that 
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there is a constant c > such that given N n = k, Vox(H(A n )\J : ) > ck. Therefore 
Vax(H(A n )) > cE [N n ], and this divided by \A n \ is bounded away from zero. Hence 
<7i > by (1.3). This completes the proof. | 

Finally we shall prove Theorem 1.4. The aim is to apply Lemma 3.3. Let the 
point processes £ t and £ (the set of locations of adsorbed particles at time t and 
at time oo, respectively) be as defined at (4.2) and (4.3). The family of variables 
(£(Q Z ), z G Z 2 ) forms a stationary random field. We need to show rapidly decaying 
correlations for this random field, and do so via Lemmas 5.4 and 5.5 below. The 
first of these (but apparently not the second) can be proved using Theorem 4.20 
from Chapter I of Liggett [7], but we take a different approach which is closer to 
that used already. 

For z G Z d , let B" z be the union of all sets B' ,y G B(z; lOp + 2). Lemmas 4.1 
and 4.3 imply an exponentially decaying tail for the diameter of the set B' , and 
hence the distribution of the diameter of B" also has an exponentially decaying tail, 
uniformly in z, i.e. 

Hmsupsupr _ MogP[B"n(Z 2 \.B(z;r)) ^ 0] < (5.6) 

Lemma 5.4 Given any finite T C Z 2 ; let jF r be the a -field generated by (£(Q Z ), z G 
T). There exist positive finite constants K',5 2 such that ifT,V are sets in Z 2 , both 
of cardinality at most 4, and the distance between them is d(T,T'), we have for all 
events F G Ty, G G Tv , 

P[F n G] - P[F]P[G] < K / exp(-5 2 d(r,r')). 

Proof. Let V be an independent copy of the Poisson process V . Suppose F G JF r 
and G G Ty>- Then F = {(£(Q z )) z£ r G R} for some Borel R C R r , and G = 
{i(Qz)z&v> G R') for some Borel R' C R r ' . 

Let Hp be the set of points of M 2 lying closer to F than to G and let Hq = ~R 2 \Hp. 
Let F* be defined like F but based on points of the Poisson process 

Q : = (V n (Hp x [0, oo))) U (V H (if G x [0, oo))). 

That is, let £ ( ^ 1 - ) be defined in the same manner as £ (eqn. (4.3)) but in terms of the 
Poisson process Q instead of V, and let F* := (Q 2 )) 2er G /?}. Similarly, let 
G* be defined like G based on points of 

Q' := (P n (H G x [0, oo))) U (V n x [0, oo))), 

that is, let ^ be defined in the same manner as £ (eqn. (4.3)) but in terms of the 
Poisson process Q! instead of V, and let G* := {(£ (2) (Q z )) 2Gr G R'}. 
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Then F* and G* are independent (since based on independent Poisson processes 
Q, Q>), and P[F*] = P[F], and P[G*\ = P[G\. Therefore 

P[F n G\ - P[F\P[G\ = P[F n G] - P[F* n G*\ 

so that 

\P[F n G] - P[F]P[G]\ < P[FAF*] + P[GAG*}. 

By Lemma 4.4, FAF* does not occur if Q y C Pf for all y G -B" and all x e T 1 . 
Likewise, GAG* does not occur if Q y C -Pf for all ?/ G -B" and all a; G T 2 . By (5.6), 
P[FAF*] and PfGAG*] both decay exponentially in d(Ti, r 2 ), uniformly over finite 
Tj C Z 2 of cardinality at most 4 and over F G JF ri , G G JFp 2 . g 

Lemma 5.5 Let JF = J^{o} fe£ jF t be the a -field generated by the variables 
£(Q z ),z G Z d \B(t). Then sup{|P[FnG] - P[F]P[G]\ : F G Fq, G G J^} decm/s 
exponentially in t. 

Proof. Let 7" be an independent copy of the Poisson process V . Suppose F G JF 
and G G JF t . Let F* be defined like F but based on points of 

(V n (B>(0; t/2) x [0, oo))) U (V n (M 2 \ ,D(0; t/2)) x [0, oo)) 

and let G* be defined like G but based on points of 

(V' D (£>(0; t/2) x [0, cx)))) U(Pn (K 2 \ D(0; t/2)) x [0, oo)). 

The precise definition of F* and G* is analogous to that used in the preceding proof. 
Then F* and G* are independent and P[F*] = P[F], and P[G*] = P[G\. Therefore, 
as in the preceding proof, 

\P[F n G] - P[F]P[G]\ < P[FAF*] + P[GAG*\. 

By Lemma 4.4, FAF* does not occur if B% C B(0; t/4). Hence by (5.6), P[FAF*] 
decays exponentially in t, uniformly over F G Tq. 

By an extension to the proof of Lemma 5.2, GAG* does not occur if 

P(0;3t/4) f)B z = Mz G A-2 n 7L d 

and the probability of this last event decays exponentially in t by (5.6). g 

Proof of Theorem 1.4- By lemmas 5.4 and 5.5, the result follows by taking 
T n = B n , in Lemma 3.3, provided we have a 2 > 0. An elementary argument shows 
that with ip z = £(Q Z ), 

|P n | _1 Var I ^ ip x J -> ^ Cov(^o,^) := 
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and the left hand side of this expression can be shown bounded away from zero by 
a similar argument to that used at the end of the proof of Theorem 1.2. Therefore 

(7 2 > 0. ■ 
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